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We discuss the problem of the equivalence between Continuous Time Random Walk (CTRW) and 
Generalized Master Equation (GME). The walker, making instantaneous jumps from one site of the 
lattice to another , resides in each site for extended times. The sojourn times have a distribution ip(t) 
that is assumed to be an inverse power law with the power index /x. We assume that the Onsager 
principle is fulfilled, and we use this assumption to establish a complete equivalence between GME 
and the Montroll- Weiss CTRW. We prove that this equivalence is confined to the case when ip(t) 
is an exponential. We argue that is so because the Montroll- Weiss CTRW, as recently proved by 
Barkai [E. Barkai, Phys. Rev. Lett. 90, 104101 (2003)], is non-stationary, thereby implying aging, 
while the Onsager principle, is valid only in the case of fully aged systems. The case of a Poissonian 
distribution of sojourn times is the only one with no aging associated to it, and consequently with 
no need to establish special initial conditions to fulfill the Onsager principle. We consider the case 
of a dichotomous fluctuation, and we prove that the Onsager principle is fulfilled for any form of 
regression to equilibrium provided that the stationary condition holds true. We set the stationary 
condition on both the CTRW and the GME, thereby creating a condition of total equivalence, 
regardless the nature of the waiting time distribution. As a consequence of this procedure we create 
a GME that it is a bona fide master equation, in spite of being non-Markovian. We note that the 
memory kernel of the GME affords information on the interaction between system of interest and 
its bath. The Poisson case yields a bath with infinitely fast fluctuations. We argue that departing 
from the Poisson form has the effect of creating a condition of infinite memory and that these results 
might be useful to shed light into the problem of how to unravel non-Markovian master equations. 
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I. INTRODUCTION 



The Onsager principle jlj is one of the basic tenets of 
statistical mechanics insofar as it establishes a connection 
between a property of equilibrium, the correlation func- 
tion of a given variable A, and the regression to equilib- 
rium of a macroscopic signal. For this reason, we judge 
the Onsager principle to be a fundamental step for the 
connection between dynamics and thermodynamics. It 
is important to stress, as clearly stated by Onsager him- 
self that this principle holds true for aged systems, 
namely systems in contact with heat reservoirs that are 
supposed to be in a condition of thermal equilibrium. If 
the regression to equilibrium is very fast, it is not so im- 
portant to ensure the equilibrium condition of the bath, 
at the moment when we begin measuring the regression 
to equilibrium of the system of interest. In fact, in the 
specific case where the bath is responsible for fluctua- 
tions that can be assumed to be white, the regression to 
equilibrium of the bath is instantaneous. The Onsager 
principle refers to a variable of interest whose dynamics 
are made stochastic by the interaction with a bath. Thus, 
when we discuss the process of regression to equilibrium, 
we have to specify if we are referring to the system of 
interest or to its bath. As earlier said, if we adopt the 
white noise approximation to describe the fluctuations 



that are responsible for the erratic motion of the variable 
of interest, then the regression to equilibrium of the reser- 
voir is virtually instantaneous and we can easily fulfill the 
condition for the validity of the Onsager principle. The 
variable of interest is characterized by a stationary corre- 
lation function $>A(ti, £2) = (A(ti)A(t2)) , which only de- 
pends on \t\ — 12 1 ■ The stochastic behavior of the variable 
of interest is caused by the interaction between system of 
interest and bath, and this kind of process is often stud- 
ied by means of the master equation method. A popular 
method to derive the master equation is the projection 
method by Zwanzig |2j|. However, this method is easy 
and convenient to use, when special initial conditions can 
be adopted, with the total distribution expressed as the 
product of the relevant, or reduced distribution, and of 
the bath distribution 3] . These special initial conditions 
have the beneficial effect of annihilating the inhomoge- 
neous term that makes the reduced equation of motion 
explicitly dependent on the initial condition If no 
initial condition of this kind is adopted, we have to con- 
sider also the inhomogeneous term and we have to wait a 
time comparable to the bath relaxation time for this term 
to disappear. This suggests that a non-Markov master 
equation with the same structure as the Zwanzig master 
equation, with no inhomogeneous term, does not nec- 
essarily refer to a stationary condition. It might refer 
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to an initial condition with system of interest and bath 
uncorrelated one from the other, a non-stationary state. 
As we shall see, to create a master equation compatible 
with the Onsager principle, in the case when the bath is 
not infinitely fast, we must establish a delicate condition 
of entanglement between system and bath. This is not 
a trivial problem, since the departure from the Poisson 
condition generates an infinitely extended memory, and 
the system-bathentanglement is the result of a rearrange- 
ment process with an infinite time scale. 

In the last few years there has been a growing interest 
for the formalism of Continuous Time Random Walk 
(CTRW) 0, because of its close connection with the 
adoption of fractional operators that, in turn, are reveal- 
ing a powerful tool to describe cooperative processes in 
condensed matter |{| . On the other hand, the connection 
between the master equation formalism and the CTRW 
walk has been discussed over the years, beginning with 
the pioneering work of Ref. [a] . According to the CTRW 
the walker travelling through an one dimensional path 
makes jumps from one site to the others after spending 
a given amount of time in the departure site. The 
distribution of sojourn times is assumed to be the same 
for all the sites, and is denoted by the function ip(t). The 
pioneering work of Ref. Q proves that the Markovian 
master equation is compatible with the CTRW, if 
is a Poissonian distribution. 

The results of the pioneering work of Ref. @ raise the 
important issue of the connection between CTRW, with 
non-Poissonian ip(t), and a non-Markovian master equa- 
tion. Apparently, this important problem is solved by 
adopting the non-Markovian master equation proposed 
by the authors of Ref. 0. We consider this theoretical 
tool to be very important, and we denote with the term 
Generalized Master Equation (GME) all transport equa- 
tions with the same structure. Recently the GME has 
been discussed by Metzler Q who argued that this equa- 
tion unifies fractional calculus and CTRW. We want to 
remark, however, that here we focus our attention on how 
to make the GME stationary, and thus compatible with 
the Onsager principle. In fact, as earlier remarked, the 
absence of an inhomogeneous term does not guarantee 
the stationary nature of the transport equation. 

We want to stress that another reason of interest of 
the GME of Kenkre, Montroll and Shlcsinger is that 
it is formally identical to the master equation that the 
authors of Ref. Q derived from a quantum mechanical 
tight-binding model, with an erratic distribution of 
energy sites. The calculation was done along the lines 
established by Zwanzig |2| using a projection operator 
similar to that used by Kenkre in an earlier work 
|10|. This means that in principle understanding the 
connection between a master equation and CTRW 
might contribute to shedding light on the intriguing 
issue of unravelling quantum mechanical master equa- 
tions. This problem becomes challenging in the case 
of non-Markovian processes. For some very recent 



references on this issue, we refer the reader to Refs. 
[Til Il2| . The purpose of this paper is not so ambitious 
as to afford direct contributions to the settlement of 
this problem. However, we focus our attention on the 
case of a two-site system that seems to be closer to the 
problem of decoherence of a q-bit than to the random 
walk over an infinite path for which the CTRW theory 
was originally designed. We address a problem that is 
closely connected to that recently discussed by Sokolov 
and Metzler [l3|. the discussion of the connection 
between the CTRW and a two-state non-markovian 
master equation. We do not discuss, however, the 
connection with the fractional calculus, an interesting 
issue addressed by the authors of Ref. [T$. Rather we 
focus on the problem of how to make both the CTRW 
and the equivalent GME compatible with the stationary 
condition and the Onsager principle. For instance, we 
define as GME the transport equation found by the 
authors of Ref. because it has the same structure 
as that proposed by Kenkre, Montroll and Shlcsinger 
7J. However, the quantum transport equation of the 
type of Ref. [!j might be non stationary, due to the fact 
it refers to an initial condition with the system totally 
disentangled from the bath. In this paper we shall show 
how to build up a genuinely stationary GME. 

The Onsager principle is a property that makes it 
possible for us to derive the correlation function from 
the GME. If we set a condition of total equivalence be- 
tween the GME and the CTRW, we find the apparently 
disconcerting result that only the Poisson statistics is 
compatible with the equivalence condition. The reason 
for this restriction is due to the fact that the departure 
from the Poisson statistics generates memory properties 
that make the GME incompatible with the Markov 
approximation. This means that the structure of the 
GME is dictated by the initial condition. If this is 
not stationary, the resulting GME is not a bona fide 
transport equation 14]. On the other hand, if the 
waiting time distribution ip(t) is not exponential, there 
are aging effects. This means that we have to leave 
the system age till it reaches the condition where the 
Onsager principle holds true. In this condition it is 
possible to establish a bona fide GME. We do it, and in 
so doing we establish a complete equivalence between 
CTRW and GME. At same time, we prove what it is 
probably the most important result of this paper. This is 
that the transition from the Poisson to the non-Poisson 
statistics has the effect of creating a memory condition 
incompatible with the Markov approximation, at least 
in principle. In fact, the Markov approximation, as we 
shall see, means that the bath is infinitely fast, and 
this condition, in turn, is proven to force the waiting 
time distribution tp(t) to be exponential. This means 
that the Markov condition is incompatible with the 
non-Poissonian nature of the waiting time distribution. 
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II. THE POISSON CASE AND THE ONSAGER 
PRINCIPLE 

In this section we prove that apparently the equiva- 
lence between the CTRW of Montroll and Weiss and the 
GME is complete only in the Poisson case. Let us express 
the GME under the following concise form: 



and 



dt 1 



I 



p(t) = - / $(t - t')Kp(t')dt', 



(1) 



where p is the m-dimensional population vector of m 
sites, K a transition matrix between the sites and <I>(i) 
the memory kernel. The prescription of the CTRW p| 
yields: 



p(t) = V / dt'Mt')*(t - t')M n p(0). (2) 

n=0 J ° 

Note that ip n (t) is the probability that n jumps occurred 
and that the last took place at time t — t' . This means 
that 



Mt) = / Ipn-iit-t')^')^', (3) 
JO 



with 



(4) 



While M is the transition matrix connecting the sites 
after one jump has occurred. The function 'S(i) is the 
probability that no jump occurs in the time interval t, 
namely, 



*(*) 



tpi(t')dt'. 



(5) 



It is evident that the time convolution structure of both 
Eq. (JTJ and Eq. J2J makes it straigthforward to derive in 
both cases the Laplace transform of p(t). Furthermore, 
Eq. @ and Eq. JSJ) makes it possible to express the 
Laplace transform of p(t), in the second case, in terms 
of the function ipi (t) only. This function is the earlier 
mentioned waiting time distribution, and from now on it 
will be denoted with the symbol ip(t). By comparing the 
Laplace transform of the GME to the Laplace transform 
of the CTRW we get 



wp(u) M - I 
1 - iJj(u) K ' 



(6) 



where I is the unity operator. We limit our discussion in 
this paper to the two-state case, where 



M = 



1 

1 



(7) 



K 



1 -1 
-1 1 



thereby turning Eq. © into 
l>(u) 



iip(u) 



(8) 



(9) 



1 - i/>(u) 

We note that the two-state master equation studied 
in this paper has to be related to the two-site CTRW. 
If we assign the value W to the right site, and the 
value — W to the left site, and we adopt a discrete 
time representation, the motion of the random walker 
corresponds to a symbolic sequence {£}, with the form 
{WWWWWWW W -W-W-W- WWWWWWW - 
W-W-W-W-W-W-W...}, which shows a signif- 
icant time persistence of both states. The waiting time 
distribution tj}(t) is the distribution of the patches filled 
with either W's or — W's. We assume a symmetric condi- 
tion. In the case when the first moment of ip(t) is finite, 
the process is stationary 0, thereby making it pos- 
sible to define the correlation function of the fluctuation 
denoted with the symbol 

It is now the right time to use the Onsager principle. 
We set the Onsager principle in the form 

Pl(t) -P2(t) 



(10) 



Pi(0)-pa(0)' 

where p\ (t) and P2 (t) are the probabilities for the random 
walker to be, at time t, in the first and second state, re- 
spectively. In his original work Onsager referred himself 
to the case of a macroscopic fluctuation that is supposed 
to regress to the vanishing equilibrium value through a 
phenomemological equation of motion. We realize this 
macroscopic fluctuation selecting a large number N of 
walkers, divided into two groups, with Ni and iV 2 walk- 
ers belonging to the first and second state, respectively. 
Then we relate this choice to the probabilities pi(t) and 
P2{t), by setting pi = Ni/N andp2 = N%/N. The regres- 
sion to equilibrium of this macroscopic fluctuation does 
not fit any phenomenological law. For this reason, we can 
refer ourselves to Eq. I|1(J|I as a generalized version of the 
Onsager principle. We plan to derive the mathematical 
expression of the law of regression to equilibrium from 
the use of this generalized version of the Onsager's re- 
gression principle, which, as shown in Section 3, is exact 
in the dichotomous case. Due to its nature, the GME can 
afford information only on p\ (t) and p2 (t) . There is no 
direct information on the correlation function $ t (t) , and 
the only possible way is through the Onsager's principle, 
which, as shown by Eq. I|1(J|) . implies that an initial off- 
equilibrium condition, pi(0) — P2(0) ^ 0, is established. 
Then, the regression to equilibrium is expected to take 
place in such a way as to establish a connection with 
the correlation function in the form of Eq. llUfl . This 
assumption, applied to Eq. JQl, yields: 



(11) 
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This equation is the non-Markovian counterpart of the 
phenomenological regression to equilibrium of the origi- 
nal work of Onsager Using the time convoluted ex- 
pression of Eq. 1|11[> and the central result of Eq. JJjJ, we 
establish a connection between the correlation function 
and the waiting time distribution tp(t), through 
their Laplace transforms, as follows: 



1 



2utp(u) 

(i-v(«)) 



(12) 



On the other hand, the correlation function is con- 
nected to ip(t ) th rough properties established by the re- 
newal theory [l5j. To properly establish this connection, 
we have to notice, with Zumofen and Klafter ^3] , that it 
is convenient to introduce another type of waiting time 
distribution that we call ip*(t). What's the connection 
between tp{t) and ip*(t)7 The waiting time distribu- 
tion ip(t) is the experimental waiting time distributon. 
It could be evaluated experimentally observing the se- 
quence {£}, and recording the time length of the laminar 
regions occupied only by W or by —W. However, we 
can imagine that a theoretical waiting time distribution 
exists, denoted by ?/>*(£), and that the sequence is 
obtained as follows. We select randomly a number t\ 
from the distribution ijj*(t). We toss a coin, and assign 
to the first laminar region, of length ti, the symbol W or 
— W, according to the coin tossing prescription. At the 
end of this laminar region, first we select, again from the 
distribution ip*(t), a number t2- This is the length of the 
second laminar region. We toss the coin again to decide 
the sign of it. It is evident that there is 50% probabil- 
ity of getting the same sign as the earlier laminar region. 
We proceed in the same way with the length and the sign 
of third laminar region, and so on. We adopt this rule 
to create the sequence {£(/)}• Thus, from the renewal 
theory [l5l | we obtain the following important result: 



{t' -t)^*{t')dt', 



(13) 



where (t) is the mean waiting time of the tp*(t)- 
distribution. It is interesting to notice that this equa- 
tion means that the second derivative of the correlation 
function is proportional to ip*(t), 



dt 2 



(t) ' 



(14) 



Zumofen and Klafter [l7j explained with clear physical 
arguments the connection between ip(t) and ip*(t). They 
prove that their Laplace transforms are related one to 
the other by 



2j){u) 

i + i){u) ' 



(15) 



Using Eq. (|14f) and Eq. (|15fl we derive a further con- 
nection between the Laplace transform of if) and the 



Laplace transform of ifi(t). In fact, by Laplace transform- 
ing Eq. I|14|) and using Eq. I|15|) , we derive a new expres- 
sion for the Laplace transform of <E»^ (i) . By equating this 
new expression to that of Eq. I|12|). we find, after some 
algebra, the following form for the Laplace transform of 

m, 



■K(o) 



+ WO) 
Using Eq. (fT^|) we prove that 



(16) 



(17) 



This equation allows us to write Eq. (|16fl under the form 



2 <*> ( M +2^y) 



(18) 



This means that ip{t) is exponential and the explicit 
forms of ip(t) and ^^(t) are given by 



and 



V>(t) = |exp(-^) 



=exp(- 7 t), 



(19) 



(20) 



respectively. 

This is the first important result of this paper. It is to 
some extent disconcerting since it seems to restrict the 
complete equivalence between GME and CTRW to the 
exponential case, with the effect of making the GME use- 
less, since in the exponential case, it turns out to be an 
ordinary, memory-less, master equation. We want to no- 
tice that Sokolov and Metzler [lj|] have recently discussed 
the derivation of a fractional transport equation totally 
equivalent to the CTRW. However, they did not discuss 
the intimately related problem of the formal equivalence 
between CTRW and GME, which is the main goal of the 
present paper. 

The second important result of this paper is the discov- 
ery of a GME that is equivalent to the CTRW, with no 
restriction to the Poisson case. This second important 
result will be illustrasted in Section 5. To derive this 
important result we have to mention a fact recently ob- 
served by Barkai 0|: the Montroll- Weiss CTRW refers 
to a non-stationary condition. This means that the On- 
sager's regression principle is invalidated. Therefore, we 
shall conclude that the disconcerting result of this section 
is due to the fact that the case of Poisson statistics is the 
only one where aging does not exist. At this stage, this 
is only a conjecture. We can notice, however, that this is 
a plausible conjecture. In fact, the exponential condition 
of Eq. p9fl when applied to Eq. @ makes <t(u) indepen- 
dent of it, thereby implying that instantaneous regression 
to the vanishing value of the memory kernel &(t). This 
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means that the fluctuations of the variable £ take place 
with its bath at equilibrium. 

We give a further support to the conjecture that ag- 
ing is the reason of the apparently disconcerting result 
of this section by referring ourselves to the work of Ref. 
. The authors of this paper found a transport equation 
with the same form as that of Eq. , with no inhomo- 
geneous term. The adoption of the Zwanzig projection 
method makes this interesting result possible, provided 
that special initial conditions are chosen. These condi- 
tions might depart from the entanglement between the 
system of interest and its bath, necessary to ensure the 
stationary condition. On the other hand, the Hamilto- 
nian treatment implies that the system exerts an action 
on its bath. Thus, the choice of an initial condition an- 
nihilating the inhomogeneous term might have the effect 
of perturbing the bath, with an ensuing long lasting re- 
equilibration process, and so a persistent non-stationary 
condition. 



III. ONSAGER'S REGRESSION FOR 
DICHOTOMOUS SIGNALS 



We note that the constraint posed in the earlier section 
on the waiting time distribution tp(t), to ensure the equiv- 
alence between the Master equation and the CTRW, has 
been based on the assumption that the Onsager principle 
is valid. The fact that the equivalence between CTRW 
and GME is restricted only to the Poissonian case, and 
so to the case of ordinary master equation, might gen- 
erate the false impression that the Onsager principle is 
valid only in the Poissonian case. It is not so. It is al- 
ready known that in the case of Gaussian statistics 
the Onsager's regression hypothesis is exact, and it holds 
true for initial excitation of any intensity. In this sec- 
tion, having in mind the two-state master equation, we 
focus our attention on the case of dichotomous statistics. 
In the dichotomous case, we show that the Onsager's re- 
gression hypothesis turns out to be exact again, for initial 
departures from equilibrium of whatsoever intensity, and 
with any type of relaxation process, with an exponential 
form, and with an inverse power law form as well, pro- 
vided that the system is ergodic. In short, in this section 
we prove that in the dichotomous case the Onsager re- 
gression principle holds true provided that the system is 
ergodic. 

The ergodic property implies that we can create a 
Gibbs ensemble in two different, but equivalent, ways. 
The first refers to the original idea of Gibbs, that the 
same system is repeated infinitely many times. This 
means that we have to generate infinitely many sequences 
using the same physical prescription. The second way is 
based on the adoption of only one sequence, denoted by 
£(f), which is supposed to be infinitely long. Then we 
define a generic trajectory with the following pre- 

scription 

£ W (t) + (21) 



In the practical case of a numerical treatment, time is 
discrete, time unity being, for instance t = 1. In this 
case the superscript s is an integer number. 

We denote by P+ and P_ the probability that the vari- 
able £ gets the value W and —W, respectively. These 
probabilities are frequencies that can be evaluated us- 
ing either the first or the second ensemble of sequences. 
We divide our ensemble of sequences into two groups, 
characterized by the the initial condition £(0) = 1, and 
£(0) = —1, respectively. We define two distinct averages 
for the two groups, denoted by (•}+ and (•)-, respectively. 
Of course, we have that: 



with 



(•)=P + (.) + + P_(.}_, 



P+ + P- = 1. 



(22) 



(23) 



Although it is not essential for the main goal of this 
paper, we make our discussion as general as possible. 
Thus, we do not assume the two probabilities to be equal. 
We do not rule out that a bias might exists, given by 



P+ - P_ = c. 



(24) 



In the ergodic, and stationary, condition that we are as- 
suming, the stationary autocorrelation function exists, 
and it is defined as 

= TTTx 77\2 = 1 3 • ( 25 ) 



l-c 2 



Expressing the total average in terms of the averages over 
the two groups, according to the prescription of Eq. J22J>, 
we make Eq. (|25ll become 



1 -c 2 
l-c 2 



(26) 



Note that on the same token the property (£) = c can be 
written as follows 



p + (0+ + P-{0- = c 



We also note that 



c + 1 



(27) 



(28) 



Using Eqs.J23 an( i (|28[) we express the correlation func- 
tion of Eq. (|26|l as follows: 



* c (t) = 



2P + (0+-c-c 2 _ (£} + ( c +l)- c ( c + l) 



l-c 2 



l-c 2 



(29) 

Finally, by dividing both numerator and denominator of 
the ratio corresponding to the last term of Eq. I|29|l by 
1 + c, we obtain 



(0- 



1 - c 



(30) 
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Eq. I|30[) can be interpreted as follows. The group of tra- 
jectories corresponding to the condition £(0) = 1 can be 
thought of as a way of creating an out of equilibrium con- 
dition, (£)+. This out of equilibrium condition undergoes 
a process of regression to equilibrium that is proportional 
to the equilibrium correlation function. Using Eq. Ij27|) . 
we can express the correlation function in terms of (£)_, 
and we reach the same conclusion. It is evident that 
the same conclusion would be reached using an arbitrary 
mixture of and (£)_, departing from the vanishing 
equilibrium value. In conclusion, the Onsager's regres- 
sion principle is fulfilled. 

We notice that the waiting time distributions that we 
are considering in this paper have an inverse power law 
nature with the power index fi. If this power index fits 
the condition \x > 2 the ergodic condition is ensured. We 
shall focus our attention on this condition, and we shall 
prove that we can make the CTRW compatible with the 
GME if we set the constraint that the system is aged 
enough as to make the Onsager principle valid. As we 
shall see, this means the adoption of a form of CTRW 
different from that of Montroll and Weiss Q , which corre- 
sponds to a condition very far from the stationary state. 
It is evident that the ergodic condition might become 
hard to fulfill in practice, with relaxation processes de- 
scribed by inverse power laws. This is in fact the case 
when aging becomes important. 



IV. AGING IN RENEWAL PROCESSES 

In Section 2 we found that the GME and the CTRW 
are equivalent only in the Poisson case. Since the use of 
the Onsager's regression assumption is the key ingredient 
used to establish this equivalence, we might be tempted 
to conclude that the Onsager principle does not hold true 
in the non-Possonian case. In Section 3, we found that it 
is not so, and that the Onsager principle is valid, provided 
that fj, > 2. Here we shall prove that the non-Poisson 
case results in aging effects, which has significant conse- 
quences also in the case /i > 2. If the CTRW used does 
not refer to the stationary (aged) condition, it cannot be 
compatible with the Onsager regression principle. This 
is the second result of this paper. The problem of ag- 
ing within the context of intermittent processes has been 
discussed recently in a very attractive paper by Barkai 
|l8j. The work of Barkai focuses on the non-stationary 
condition, fj, < 2, where no invariant measure exists, and 
aging lasts forever. In this case the Onsager regression 
principle is violated for reasons explicitly pointed out by 
Onsager himself. In fact, Onsager repeatedly mentioned 
that his principle has to be applied to aged systems. 
The condition of aged systems implies a state of ther- 
modynamic equilibrium. In the case fj, < 2 equilibrium 
is never reached, not even asymptotically. Thus, in this 
case the Onsager's principle is always violated. In Sec- 
tion 4A, using the renewal theory [16|. we shall study the 
case fj, > 2, and we shall discuss two limiting conditions, 



corresponding to the birth and the death of the system, 
respectively. It is interesting to remark that death, mean- 
ing an ordinary thermodynamical condition, takes place 
eventually at the end of an infinitely long aging process. 

A. Young and aged systems 

First of all, we should define, as Barkai does |T§|. the 
waiting time distribution tpt a (t)- This means that we 
establish non-stationary conditions at time t a < 0, and 
we begin our observation at time t = 0. The waiting time 
distribution depends on t a . The naive conviction that 
the waiting time distribution is given by tp(t), actually 
is based on the assumption that t a = 0. This means, 
in fact, that we are considering a set of random walkers 
and that at time t = all of them begin their sojourn in 
the laminar region. This leads immediately to tpt a =o(t) — 
tp(t). It is straighforward to evaluate also the distribution 
of sojourn times corresponding to t a = — oo. We denote 
this distribution with the symbol tpoo(t)- 

This is the stationary case, corresponding to the follow- 
ing procedure. In the stationary condition, the probabil- 
ity of selecting a laminar zone of length T, by a random 
choice, is rxip(T)dT ', where (r) is the mean length of a 
laminar zone. The probability density of observing the 
first change of laminar phase after a time t, being in a 
laminar zone of length T, is Q(T — t)^. Consequently for 
the probability density of having the first change of lam- 
inar phase at time t (i.e. ipoo(t) ), we have, integrating 
over all possible T; 

oo 

iM*) = T^J dTT^(T)^0(T-t) 



oo 

= ±-J dT^(T). (31) 
t 

Note that due to the renewal theory, the distribution 
V'oo (i) concerns only the time that we have to wait to de- 
tect the first event. After the first event we have a total 
rejuvination. In fact, measuring the time at the moment 
of the first jump is equivalent to beginning the measure- 
ment process at the precise moment when the walker en- 
ters the laminar region. Thus the Laplace transform of 
p(t) is given by 

,* , s "0oo(w)^(u)M N . . . , 
p u) = » / + ™ K > \ p . 32 
1 - 4>(u)M 

This result is derived from Eq. (J2J) as follows. First of all, 
we evaluate the Laplace transform of Eq. @ ■ We obtain 

p(u) = (#(7i) + *(m)v3(u)M + *(w)^ 2 M 2 + ...)p(0) 

= f>M)"p(0) = y. p(o). 
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Then, we replace the probability of occurrence, or of non 
occurrence, of the first event, called %p(t) and ^(t), re- 
spectively, with the corresponding aged quantities. These 
are called tpoo(t) and ^^(t), respectively. The function 
V>oo(i) is given by Eq. l(3*Tl) and the function ^oo{t) is 
given by 



4>oo(t')dt'. 



(33) 



All this yields Eq. 

We expect that the inverse Laplace transform of Eq. 
|J3j5}, is a function of t that asymptotically will become 
equivalent to the asymptotic value of the inverse Laplace 
transform of Eq. I|32(l . What about the Onsager prescrip- 
tion of Eq. (|10|l ? If we adopt the GME corresponding to 
the young condition, the adoption of Eq. I|10() would be 
equivalent, as far as the left-hand-side of this equation 
is concerned, to adopting the non-stationary condition. 
We have to remark, in fact, that the GME corresponding 
to the Montroll- Weiss CTRW is not a bona fide mas- 
ter equation. The importance of a bona fide transport 
equation has been pointed out in a remarkable paper by 
Fox 0| . A master equation is a bona fide master equa- 
tion when it can be used with any initial condition (like 
the bona fide Fokker-Planck equation of Ref. ^4| ) • The 
master equation of Section 3, on the contrary implies 
the choice of only one initial condition. The next subsec- 
tion is devoted to sheding more light into this important 
aspect. 



B. A theory for systems of any age 

The purpose of this subsection is to shed further light 
into the aging problem, and into the ensuing conflict with 
the Onsager postulate. At the same time we shall derive 
an expression for the distribution of the first exit times 
valid for any age, and not only for t a = and t a = — oo, 
which denote observation taking place at the moment 
of birth and death of the dynamic process under study, 
respectively. This discussion is based on the following 
dynamic model. A variable x moves in the interval / = 
[0, 1] according to the following prescription 



(34) 



When it reaches the point x = 1 it is injected back with 
uniform distribution, thereby producing another laminar 
region. We assign alternated signs to the sequel of lami- 
nar regions. The sojourn in a laminar region with a given 
sign is equivalent to sojourn in one of the two states dis- 
cussed in Section 2. The resulting waiting time distri- 
bution coincides with the one earlier called ip(t). The 
explicit expression for ip(t) is obtained by expressing the 
exit time f as a function of the initial condition xq . Then 
we have to assume 



The choice of a uniform of back injection process implies 
p(x ) = 1. This is the condition behind the CTRW of 
Montroll and Weiss. All this results into 



1&(t) = (/i-l) 



(T + t)v ' 



with 



and 



(z-1) 



T = 



1 



X(z-l) 



(36) 



(37) 



(38) 



We propose for the aging process a calculation proce- 
dure different from that adopted by Barkai 0]. As we 
shall see, this procedure yields the same results as those 
of Barkai. As done in Ref. [2(J, we discuss the equation of 
motion for the probability density p(x,t). This is given 
by 

-p( x ,t) = - — (\x*p(x,t))+C(t). (39) 

The time evolution of the distribution function is de- 
scribed by the following formula 



p{x,t) 



C(T + t a ) 



_ t [1 + l/a(t - t)^ 1 /"]^" 



G?T, 



where we have used the following new parameter: 



(z-1) 



(40) 



(41) 



We have assumed that the flat initial distribution is 
assigned at time t = —t a . We have neglected the contri- 
bution [1 + l/a(t + t a )( 1 / a )]~'-- 1+a \ necessary to recover 
at t = —t a , the flat initial distribution. On the basis of 
the approach that we shall detail hereby to derive the 
distribution of the first exit times, it is straightforward 
to check that this term yields negligible contributions. 
Using the property 



p(x, t — 0)dx — ipt a (t)dt, 



(42) 



where p(x, t = 0) is given by Eq. l|40|l . with t = 0, we can 
evaluate the distribution function of the first exit times 
in general, and not only for t a — and t a = — oo. It is 
interesting to remark that the infinitely old distribution 
V'oo(^) of Section 4 A, is easily obtained by noticing that 



p{x, oo) = (2 — z)jx 



2-1 



(43) 



p(xo)dxo = ijj(t)dt. 



(35) 



Using Eas. l|43|l and (|42(l . after some algebra, we rederive 
the distribution of Eq. I|31|) . 

We can also establish a connection with the case dis- 
cussed by Barkai [l8j, namely the case z > 2. In the case 
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of z > 2 =>■ a < 1, the long time behavior of the function 
C(t) is 



(44) 



We now use the prescription of Eq. I|42|) . A close inspec- 
tion of Eqs. (|40(1 and (|44Jl reveals that it is not possible 
to give a close formula for the function tp ta (t) . Then, we 
must study the two cases t a 3> t and t a <C t, separately. 



• t a > t: 

p{x,t = Q) 



t a ) 



-(1- 



[1 - l/arx 1 l a ] 1+a 



dr 



= t: 



(1 + 2/) 



-(l-c 



.! [i/i a -y/aa; 1/Q ] 1+0 

1 (i-y)-a-") 

[l/t a + y/as 1 /*]!-** 
1 



dy 



1 



o [l/t„ + y/oarVaJi+a 



ta J [l + y/at a x l / a Y+ a 

.-(l-a) 



,1/a 



(45) 



In Eq. we have carried out the operation lim ta ^ 00 , 
while keeping x ^ fixed. The next step is based on 
the use of Eq. g2J>, of x = (1 + i/a) -0 and of | = 
(1 + t/ay {a+1) . This allows us to obtain 

^*a(*)lt »t 



'•a 



• t a < f: 

p{x,t 



0) 



to [1 - l/ara; 1 /"]^" 

3 (T + t„)~ (1_o) dT 



Thus, in this condition we have 

±a 

< n 



(47) 



t l+a 



Let us note that Eqs. I|4t)|) and (I48|) correctly reproduce 
the behavior predicted by Ref. 18] ) in the corresponding 
limits by means of the formula: 



sin(a7r) t" 



t<*(t + t a )' 



(49) 



where a = —^j . In conclusion, our way of proceeding fits 
the conclusions of Barkai [l^. As earlier pointed out, 
in the condition \x < 2 the Onsager principle is always 
violated. This is expected, since it is known (2l| that the 
non-stationary condition /i < 2 is incompatible with the 
existence of the stationary correlation function $f (i). 



V. THE STATIONARY MASTER EQUATION 

In this section we prove that the GME can be made 
equivalent to the CTRW, with no restriction to the Pois- 
son statistics, if we set the stationary condition in both of 
them. The deep difference between stationary and non- 
stationary condition has been discussed years ago by Zu- 
mofen and Klafter in two pioneering papers |23l . This 
means that we have to compare the GME of Eq. to 
the stationary CTRW, which yields 



p M = (f M M/ /r ( f ( f )p(«). (so) 



l-i>(u)M 

where, according to the stationary prescription |2fij : 

1 f°° 

iM*) = T^y J *W)> ( 51 ) 

as also explained in Section 4A, and (r) is used again to 
denote the mean waiting time of the distribution tp(t). 
For M we make the same choice as in Eq. {7J). For cal- 
culation convenience, we adopt for K the more general 
symmetric form: 



K 



x y 



(52) 



(46) From (JTJ we get: 



p(u) 



1 



ul + §(u)K 



p(0). 



(53) 



By equating (|53|l and l|5U|l we derive 



(48) ul + <5>{u)K 



+ ^oo(w)*(ii)M J2 i>{u) n M n . 

(54) 

By multiplying both sides by ul + $(u)K and using the 
properties M 2n = I and M 2n+1 = M, we obtain: 



/ =*„(«)(«/ + $(u)K) + Voo(u)*(w) ^—M 

\1 — tp(u) 2 

' mk i ^ u)u I I k] 



1 - V(w) 2 



l-V(w) 2 l-V(u) 2 



(55) 



With some algebra we extract from Eq. (|55|l the following 
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set of equations for x, y and z: 



V»oo(«)*(«) ^ 

=*oo (u)$(u)y+ 

^oo(u)*(u) f - 



-z + 



=* 00 («)$(«)y+ 



-a: + 



1 + $(u)z)+ 



lf)(u)u 


^(u)d(u) 


1 -*0(u) 2 


T \ J 




1 $(«)*(«) 


1 -^(li) 2 


1 - VXw) 2 


U 


! $(«)*(«) 


1 - ^(w) 2 


1 -i/>(u) 2 


%p(u)u 


$(«)*(«) 


l - v3(w) 2 


1 -i/)(u) 2 



The solution of this set of equations is given by x — z, y = 
—x and 



it(l — ^(u)) 



$(w)(-2 + u(t) + (2 + ^(t))^)) 



(56) 



So, the GME and the stationary CTRW are compatible 
if K has the same form as that of Eq. (JSJ and 



u(l — ip( u )) 



(-2 + u(t) + (2 + u(t»V(u)) 



(57) 



Using the Onsager's principle in the form of Eq. (|10|l and 
consequently the time convoluted equation of Eq. 
we obtain 



*£(«) 



1 



It + 2$(u) 



(58) 



Finally, using for <t(u) the expression of Eq. I|57[) and 
for $^(m) the expression afforded by the renewal theory, 
we get an identity in ip(u) (the Laplace transform of the 
waiting times distribution). This means that the aged 
CTRW is totally equivalent to the GME with the memory 
kernel given by Eq. Q57[l. 



VI. CONCLUDING REMARKS 

This paper affords a significant contribution to under- 
standing the role of aging in systems departing from Pois- 
son statistics. The Onsager regression principle rests on 
the second principle of thermodynamics, and it applies to 
aged systems. However, non-Poissonian statistics might 
imply a very slow regression to equilibrium, so as to cre- 
ate conditions, correctly depicted by the Montroll- Weiss 



CTRW, but incompatible with the Onsager regression 
postulate. This might be a condition different from the 
out of equilibrium thermodynamics, and represent in- 
stead a state of matter intermediate between the dynamic 
and the thermodynamic condition. This new condition 
has been recently proposed to explain the emergency of 
life on earth and has been denoted with the appealing 
term of living state of matter [22| . 

The master equation of Section 5 corresponds to the 
aged systems where the Onsager principle holds true. 
The Montroll- Weiss CTRW and the corresponding mas- 
ter equation refers to the dynamics of a young system. 
From a technical point of view one might wonder if it 
is possible to produce a master equation, or a CTRW, 
reflecting correctly the process of aging. The young and 
the old sequence reflect two different initial conditions of 
the system. Moving from the young condition the system 
will eventually relax to the equilibrium, but this relax- 
ation will not follow the correlation function behavior. 

Eventually, the system will settle at equilibrium. 
However, to establish, at least ideally, the validity of 
the Onsager regression principle, we have to set some 
suitable out of equilibrium initial conditions. If we do 
select initial conditions that do not fit the stationary 
prescription |23l |24j, the correlation function ^^(t) 
cannot be defined, and consequently we are forced to de- 
termine the memory kernel of the GME of Eq. Q using 
the non-stationary prescription of Eq. 0. Ncvcrthless, 
this equation will produce an exact time evolution. Let 
p(0) evolve in time till it reaches the new condition 
p(T). If we adopt this as a new initial condition, namely 
we set the new origin of time at t — T, and we apply 
again Eq. to predict the ensuing time evolution, we 
get a wrong prediction. In fact, the GME equivalent 
to the CTRW of Montroll and Weiss is not a bona fide 
master equation 0]. We should build up a new GME 
equation with the distribution of the first jumps given 
by ipt a =-T{t)- If, on the contrary, we select an initial out 
of equilibrium condition, as done in Section 3, which is 
compatible with the bath being at equilibrium, then the 
Onsager principle is fulfilled, the GME with the kernel 
of Eq. (|57|l is a bona fide master equation and can be 
used with any kind of initial conditions concerning p. 

To derive the results here illustrated, we significantly 
benefitted from the important work of Refs. 0, Il8| . 
However, we want to point out that the material deriva- 
tive of Ref. 0, although more advantageous than the 
GME to study the influence of an external force field, is 
equivalent to the GME of Ref. 0, with the condition of 
Eq. © rather than the stationary condition of Eq. I|57|) . 
and, consequently, as proved in this paper, it violates the 
Onsager principle. As to the work of Barkai on aging 
[ls|. we would like to note that it is confined to the 
case fj, < 2, a condition incompatible with the stationary 
master equation found in this paper, namely the master 
equation with the memory kernel of Eq. I|57|) . Barkai 
reach the conclusion that the aging effect for [i > 2 are 
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"confined to a certain time window" . We think that this 
paper has the effect of suggesting the plausible conjecture 
that actually this time window is of infinite size. Let us 
see why. We notice that in a sense the memory kernel of 
Eq. ifSTI) affords a way to measure the aging process. In 
fact, if the waiting time distribution ip(t) is exponential, 
the memory kernel is equivalent to that produced by a 
bath that would regress to equilibrium instantanesously, 
with no aging effect. The aging effect, from within our 
perspective, can also be interpreted as a regression to 
equilibrium of a system of interest, under the influence 
of a bath that is regressing, too, to equilibrium. In the 
Poisson case, the bath regression to equilibrium would 
be instantaneous, thereby yielding no aging. In the case 
of ji < oo we have a memory lasting forever. This is 
evident when /i < 3, since in this case the adoption 
of Eq. (|57|l proves that the Laplace transform of the 
memory kernel, tends to vanish for u — > 0. This 

is an indication of infinitely persistent oscillations that 
are clearly incompatible with the Markov approximation. 

Apparently, the condition fj, > 3 seems to be compat- 
ible with the Markov approximation, given the fact that 
in that case, tends to a finite value for u tending to 
0. We note that, while the analytical expression of $(£) is 
at the moment not yet known to us, the analytical form 
of $^ (t) is under our total control, and it is a an inverse 
power law with index f3 = fj, — 2. The free diffusion pro- 
cess resulting from the fluctuations of £ (t) is known. The 
pioneering work of Ref. 17] proves that even in the case 
/i > 3 the Gaussian condition is not exactly realized. The 
central part of the diffusion pdf is the ordinary Gaussian 
diffusion, but at large distances slow tails with an inverse 
power law appear. These tails become weaker and weaker 
upon increase of fj, and might be annihilated by even weak 
external fluctuations, with the ensuing annihilation of ag- 
ing. However, there are deep physical motivations to fully 
understand the theoretical issues generated by the aging 
process. A relevant example of this theoretical request is 
given by the physics of glassy systems, which are charac- 
terized by aging and consequently by a violation of the 
ordinary fluctuation-dissipation condition, which is inti- 
mately related to the validity of the Onsager principle 
|26j. In other words, we have to do further research work 
to understand how robust the non-Poissonian nature of 
t/j(t) is, and to establish the intensity of environemental 
fluctuations that are expected to produce a truncation of 
the inverse power law nature of ip(t). This can be done 
along the lines of Ref. |2]J. This research work can be 



carried at a much deeper level, as a consequence of the 
fundamental result of Eq. lj57|) . In fact, this equation de- 
fines the memory kernel of the GME, in terms of the 
waiting function ip(t). We assume that the function ip(t) 
can be derived from the experimental observation as in 
the case of the blinking quantum dots of Ref. (2{|. Thus, 
with the help of Eq. (|57|l we can draw from the experi- 
mental function ip(t) information on the bath responsible 
for the intermittent character of the process under study. 
This information will help the foundation of a dynamic 
model, and consequently will make possible to discuss 
the robustness against noise of the intermittent process 
under study, from a perspective more realistic than in 
the earlier work, for instance, the work of Ref. |27|. 

Let us summarize the main results of this paper. The 
functions &(t) and ip(t) are the- main mathematical prop- 
erties behind the GME and CTRW, respectively, the for- 
mer affording indications on the complex dynamics of 
the irrelevant variables and the latter illustrating the na- 
ture of the intermittent process that can be observed as 
a result of a special experimental detection |2f| . The de- 
viation of ^(t) from the exponential condition is equiv- 
alent to forcing Q(t) depart dramatically from the white 
noise condition. The non-Poisson nature of ip(t) gener- 
ates memory infinitely extended in time, and aging with 
no finite time scale as well, no matter how big the scal- 
ing parameter /x, provided that /i < oo. In fact, if a 
finite time scale for both aging and memory existed, at a 
much larger time scale, the function would be equiv- 
alent to a delta of Dirac, the bath fluctuations would be 
white noise, the bath would be at equilibrium, and, of 
course, tp(t) would be exponential. This is equivalent to 
proving that the non-exponential nature of ip(t) yields 
aging effects that are not limited to the non-stationary 
case /x < 2, even if in this case these aging effects seem 
to be more natural. The condition \x > 3 is compatible 
with the stationary condition. However, a bath gener- 
ating this condition, through the key result of Eq. I|57|) . 
is characterized by dynamic properties with memory in- 
finitely extended in time. 
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